Splitting: Tanaka's SDE revisited 
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What follows is my attempt to understand a set of ideas being developed by Boris 
'q^'. Tsirelson. I do this by studying a specific, and I hope interesting, example. 
Q>^ ' Tanaka's SDE is one of the easiest examples of a stochastic differential equation 

0\ • with no strong solution. Suppose (X^; t > O) is a real- valued Brownian motion starting 

^ • from zero and we put Bt = sgn{X^)dXs then B is also a Brownian motion and 

Q . Tanaka's SDE 

(1) Xt= I sgn{Xs)dBs, 



> 
in 



is satisfied. But the trajectory of B does not determine that of X. Recall Tanaka's 
formula 

(2) |Xi|= / sgn{X,)dX, + Lt, 

Jo 

where (L^; t > O) is the local time process of X at zero. We find 

(3) \Xt\ =Bt + sup(-E,) 

s<t 

but B does not tell us the signs of the excursions from zero made by X. 

In a discrete time framework things, work out differently. Let > O) be the 

. symmetric nearest neighbour random walk on the integers. Define sgn{a) to be +1 if 
"Th ! a > and — 1 if a < 0, and then let Zq = and Z^+i — Zn = sgn{Xn)[Xn+i — Xn) 
I then n > 0) is again a symmetric random walk and we may write a discrete 
version of equation (|^): 

>< : (4) ^n = Yl ^aniXk) {Zk+i - Zu) . 

: 

The equations @) and (H) have discrete time versions: 

n-l 

(5) \Xn + \\-\ = ^sgn{Xk){Xk+i - Xk) + U 



k=0 



where Lq = and for n > 1 we define L„ = Ylk=Q ^( ^^^i 
(6) |X„ + i| - i = Z„ + sup(-Zfc). 

k<n 

The halves appear because of the lack of symmetry in our definition of sgn- it is 
not something to worry about. This time Z does determine X: the information 
about whether Xn is below or above — ^- which at first sight appears to be missing 
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in (^)- is coded in the following way. Find the last time r G {0, 1, . . . ,n} such that 
Zr = — sup;j<^(— Zfc). If this r is even then Xn + ^ is positive while if it is odd then 
Xn + ^ is negative. 

This note is concerned with understanding what happens to this precious informa- 
tion about the sign of X when we try to obtain the continuous-time model by taking 
scaling limits of the discrete model. This is inspired by work of Boris Tsirelson on the 
spectra of noises and stability- see 0,0, and 0. 

One way to understand why the information about the signs does not survive the 
passage to the limit is to observe that it is noise sensitive. Instead of one copy of 
the random walk take a pair {Z,Z') that are p-correlated where p G (0,1). This 
means that together they form a nearest-neigbour random walk on the lattice with 
E[(Zfc+i — Zk){Z'^j^^ — Zl)~\ = p. We think of Z' as being a perturbation of Z. Now 
construct X and X' from Z and Z' so that equation and its prime version hold. As 
n becomes large ( as it will when we try to take scaling limits) we find that sgn{Xn) 
and sgn{X'jJ become uncorrelated no matter how strong the (fixed) correlation p 
between Z and Z' is. Thus, in a certain sense, sgn{Xn) is asymptotically sensitive to 
small perturbations of Z. 

The discussion of the previous paragraph, although very elementary, is (a variant 
of) the observation that eventually led Tsirelson to profound results in the theory of 
filtrations . There good account of the story in . 

Next instead of constructing Z' by perturbing the whole path of Z we may only 
perturb some sections. More precisely let A C [0, 1] be a finite union of closed intervals 
with (to be safe) dyadic rational end-points. Fix p G (0, 1). For each n construct a 
random walk ((Z^, Z^); < k < n) on with 



E[(Zfc+i - Zfc)(Z^_^i - Zl 



A;2-" G A 
/t2-" G A". 



Now as before construct X and X' and consider the correlation of sgn{Xn) and 
sgn{X'^). This time as n tends towards infinity we obtain a nontrivial limit which we 
denote by 0(p, A). 

We write the Wiener chaos expansion of any random variable belonging to C^{B) 
in the form 

fo+ f h{s)dBs + [ r Usi, S2)dBs,dBs, + ... 



JQ JO 



Then we construct a finite measure on U„>o{(si, S2, ■ ■ ■ , s„) G [0, l]"'\si < S2 ■ ■ ■ < s„} 
having density |/„(si, ... , s„)p with respect to Lebesgue measure. We call this the 
spectral measure of the random variable whose chaos expansion we used. If we start 
with variable having £^-norm equal to one this measure is a probability measure- and 
we can think of it as determining the law of a finite random subset S of [0, 1]. Thus 
there is a probability |/oP that S is empty, a probability \ fi{s)\'^ds that it contains a 
single point lying in (s, s + ds) and so on. Suppose for a moment that the information 
on signs did survive in the limit, and that the X satisfying Tanaka's SDE was some 
functional of the Brownian motion B. Then we could apply this construction to 
sgn{Xi)- and obtain a random subset S. If A is a fixed subset of [0, 1] once more 



then let \S (1 A\ denote the number of points of S lying in A. Then it is reasonable to 
expect that A) would be given by 

(7) A)=E [pl^^^l] . 

The E appearing here is with respect to the law of S which does not live on the same 
probability space as X and B. 

The discussion of the preceding paragraph is based on a false premise, but nev- 
ertheless there is a random subset- still denoted by S- such that equation (|^ holds. 
This subset possess, with probability one, an infinite number of elements. If we take 
A = [0,1] then 15* fl is infinity and p°° = (by definition if you like!). Thus 
(f){p, [0, 1]) is for any p- this is just the sensitivity to noise property with which we 
began. In what follows we examine the law of this S more closely. Not surprisingly 
the Wiener chaos expansion is our principle tool. 

For any x > and t G (0, 1) let m(^t,x) denote the spectral measure of 

sup {Bt - Bh) < X 
/ie[t,i] 

Note that the total mass of this measure is just P (sup^^^jj — B^) < x) < 1, 
but we will nevertheless speak of a random set S having distribution m(^t,x)- This 
subset is supported on [t, 1]. Let qh{x, dy) denote the (defective) transition probability 
distributions of Brownian motion killed on hitting 0. 

Lemma 1. Suppose that < s < t < 1 and that S is distributed according to m(^s,x)- 
Then the subset S fl [t, 1] is distributed according to 

qt-s{x,dy)m(^t,y)- 

Now suppose that {Xt{dx);0 < t < l) is an entrance law for killed Brownian 
motion; thus, for any < s < t < 1, 



Mdy) = J qt^six,dy)Xs{dx). 
We may define a family of measures ^(i,A) for t G (0, 1) via 

m(t,x) = / Xt{dy)mi^t,y), 



and by virtue of the lemma they have a certain consistency property- that is - if 5* is 
distributed according to ?71(s,a) then for any t > s the intersection Sn[t, 1] is distributed 
according to m{t,\)- Note that the total mass of m(t,\) is the 'probability' that the 
killed BM survives to time 1 when it is started according to A- this does not depend 
on t. Because of this consistency property there is a random set S C (0, 1] whose 
distribution we denote by m\ whose intersection with [t, 1] has distribution m{t^\) for 
any t. This S may be infinite- there is the possibility of being an accumulation 
point. 



From this point on we will take (A^; < t < l) to be a multiple of the entrance law 
for the Ito excursion measure for the positive excursions of Brownian motion. Choose 
this multiple so that mx becomes a probability measure. More explicitly we have 



(8) 



y>0. 



Recall that a random variable is said to be arc-sine distributed if it has distribution 

dt 



s{dt) = 



:l[0,l](i) dt. 



7r^t{l-t) 

The time at which a BM attains its minimum between times an 1 is so distributed. 

Theorem 2. The limits (f){p,A) exist and admit the following description. Take two 
random subsets Si and S2 distributed according to mx and a [0,l]-valued random 
variable g with the arc-sine distribution . Suppose that Si, S2 and g are independent. 
Take 

S = g{l-Si)u{{l-g)S2 + g), 

then for all p and A 

<t>{p,A)^¥.[p\'^^\]. 
Proof of Lemma. Begin by writing 

1 I sup (S, - S^) < X ) = 1 I sup {Bs -Bh)<x 
\he[s,\] J \he[s,t] 

X 1 I sup {Bt -Bh)<x + Bt-Bs 
\he[t,i] 

Condition on (Br;r < t) and then replace the second factor with its Wiener chaos 
expansion and so obtain an expansion of which the typical term is 



/I f'hi rhk-i 

dBh, dBh, ... 



dB, 



1 sup {Bs - Bh) <x \ fk{t, x + Bt- Bs\hi, ... ,hk). 

\hels,t] J 

We now replace each integrand by its chaos expansion- this must simply result in the 
chaos expansion of 



1 sup (B, -Bh) <x 
\hels,i] 

On comparing the two expansions it may be seen that if S is distributed according 
to m(s,x) then S fl [t, 1] contains exactly k points at positions {hi, hi + dhi) through 
{hk, hk + dhk) with probability 



E 



1 I sup {Bs - Bh) <x] \ fk{t, x + Bt- Bs\hi, 
he[s,t] 



hk 



dhi . . . dhk, 



but since \fk{t, y\hi, . . . , hk)\'^dhi . . . dhk is just the probability distribution of S under 
m(t,y) we are done. □ 

Proof of Theorem. Stage 1. Fix an admissible subset A. For each n consider the 
correlated random walk ((Zfe,Zj[.);0 <k<ri). There is the usual weak convergence 
in the space of continuous R^-valued paths to a process ((St,SQ;0 <t<l), each 
component of which forms a one-dimensional Brownian motion and their co- variation 
is simply: 



dBtdB'^ 



pdt t e A 
dt t e A". 



Let g be the time at which B attains its minimum between times and 1, and similarly 
define g'. Now the correlation of sgn{Xn) and sgn{X[^) can be split into the sum of 
two contributions. One arises when the random walks Z and Z' attain their minimum 
(between times and n ) values simultaneously - in this case sgn^X^) and sgn{X'^ 
are equal. The remaining contribution tends to zero for large n and so the limits 
A) exist and are given by 

Stage 2. The two random times g and g' can only be equal if their common value lies 
in one of the components of A'^. For each such component we consider the common 
time at which B and B' attain their respective minimum (over that component) and 
compute the probability that this is actually the global minimum of both Brownian 
motions. We obtain 

/» /*oo /*oo 

E [1(9=9')]=/ / / p(«t,^^t;rft,c?2/i,%2)m(„,,y,)[pl^^^l]m(i_„,,y,)[pl(i-^)^^l], 
Ja<^ Jo Jo 



where 



Uf — sup{h < t : h E A}, 
vt = ini{h >t:he A}, 



and p{u, v; dt, dyi, dy2) is the law of the triple 



{g{u,v),Bu - I{u,v),By - I{u,v)), 

g{u,v) denoting the time at which B attains its minimum I{u,v) — mi{h e [u,v] : 
Bh]. 

Stage 3. By virtue of the scaling properties of BM we have 



m 



. . [ol^nAn _ r Ut+{l_t)S)nAU 



l -777/ [ J*^^-^^^ 1 

m{l-u,y)\_P \ ~ ^{{t-u)/t,y/y^)VP J" 



A well-known exercise (Revuz and Yor chapter XII) confirms that 



p{u,v -,(11,(1111,(1112) = —\t~u{dyi)\v~t{dy2) 

71 



= s{dt)\t-u)/tidyi/ Vt)\^_t)/{i-t)idy2/ - t). 

Putting these into the formula obtained in the previous section and recalling the 
definition of we obtain the desired result: 

J A'' 

□ 

It is possible to generalise the family of measures m(^t,x) from which we obtained 
rrix- Starting from a bounded function / defined on M+ we may expand 

f{Bi -Bt + x)l{ sup {Bt -Bh)<x 
y/ie[t,i] 

and whence construct a measure m{. I would like to know when such measures 
corresponding to different / are equivalent and in this case how to compute the Radon- 
Nikodym density. This is part of the problem of obtaining the spectral resolution (see 
1^) of the noise of splitting. This is a noise richer than white noise: in addition to the 
increments of a Brownian motion B it carries a countable collection of independent 
Bernoulli random variables which are attached to the local minima of B. 
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